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Unidimensional model of the ad-atom diffusion on a substrate submitted to a standing 
acoustic wave II. Solutions of the ad-atom motion equation 

C. Taillan,i'2 N. Conibc/-^ and J. Morilloi'^ 

'^CNRS; CEMES (Centre dElaboration des Matenaux et d' Etudes 

Structurales); BP 94347, 29 rue J. Marvig, F-31055 Toulouse, France. 

'^ Universite de Toulouse; UPS; F-31055 Toulouse, France 

The ad-atom dynamic equation, a Langevin type equation is analyzed and solved using some 
non- linear analytical and numerical tools. We noticeably show that the effect of the surface acoustic 
wave is to induce an effective potential that governs the diffusion of the ad-atom: the minima of this 
effective potential correspond to the preferential sites in which the ad-atom spends more time. The 
strength of this effective potential is compared to the destructuring role of the thermal diffusion and 
to the crystalline potential induced by the substrate. 



I. INTRODUCTION 

The self-organization of materials at the nanoscale 
is a promising way to avoid the expensive lithogra- 
phy stage in the conception of the semi-conductors de- 
vices. Common self-assembling techniques rely on the 
Stranski-Krastanov growth modeSfii^ buried dislocations 
networks)^ or nano-patterned substrates i^"— An alter- 
native approach to self- assemble materials at the nano- 
scale, the dynamic substrate structuring effect has been 
recently proposedi^ In this approach, a Standing Acous- 
tic Wave (StAW) governs the diffusion of the ad-atoms 
on a substrate. In the first paper of this serieSfiS we have 
established an unidimensional model of an ad-atom dif- 
fusing on a substrate submitted to a StAW and derived 
the ad-atom motion equation, a Generalized Langevin 
equation: 



(Px 



dx 



dt^ +/ l{<t),x{t'),t-t') — {t')dt' = 

~^^{x)+m+FsAw{x,t). (1) 



The first left hand side (l.h.s.) term of Eq. ([T|) is the 
usual incrtial term. The second l.h.s. is a retarded fric- 
tion force with memory kernel ^{x{t), x{t'), t — t') where 
to corresponds to the time when the StAW production 
mechanism is switched on: since j{x{t),x{t'),t — t') is a 
decaying function of |f — i'|, we can fix ig £^t — oo in the 
integral without loss of generality. The first right hand 
side (r.h.s.) term is an effective inter-atomic (substrate- 
adatom) periodic potential force derived from the effec- 
tive potential $eff(a;), the second r.h.s. term ^{t) is a 
stochastic force and the last one FsAw{x,t) is an effec- 
tive force induced by the StAW.— 

The goal of this paper is to study the solutions of 
Eq. ^ and to evidence the structuring role of the StAW 
on the adatom diffusion, through the effective Fsaw{x, t) 
force. In Ref. [lO we showed that FsAw{x,t) reads 
FsAW cos{kx -\- (p) cos{u;t) with k and to the wave vec- 
tor and angular frequency of the StAW. However, due 
to the precise nature of the ad-atom-substrate interac- 
tions, FsAwix^t) experiences also some variation at the 



substrate lattice parameter scale: the proportionality fac- 
tor FsAW and the phase factor ip vary as a function of 
the exact position of the ad-atom in between two suc- 
cessive atomic substrate potential wells. Because we fo- 
cus here on the structuring role of the StAW, we will 
only consider its large scale variation, i.e. FsAwix,t) = 
FsAW sin(fca;) cos(wt) with a constant Fsaw value and 
a constant phase factor that we fix at — 7r/2 for con- 
venience, eluding thus the possibility for this force to 
vary on the substrate lattice parameter scale: such de- 
pendence essentially affects the dynamics of the ad-atom 
and weakly the structuring effects of the StAW»^ The 
detailed study of the dynamics of the ad-atom will be 
reported elsewhere. 

Concerning the stochastic force ^(i) and the memory 
kernel 'y{x, x' , t — t'), which are related by the fluctuation- 
dissipation theoremjiSiii 



{m^{t + T))=kBTj{x{t),x{t + T), 



(2) 



with ks the Boltzmann constant and T the temperature, 
their properties are very interaction-potential model de- 
pendant and thus cannot be studied on a general ground. 
Though we have established their analytical expressions 
in a very peculiar potential case in the preceeding pa- 
per^i^ we will here use the more standard and general 
model of a centered gaussian noise for the stochastic 
force, with correlation time Tc J^"— 

• The stochastic force is then fully characterized by 
its zero mean value and autocorrelation function: 



(eW) = o, 



{amt+r)) = D- 



-kl/r. 



and the memory kernel is given by (Eq. (O): 

7(a;(t),x(i + T),r) = — e -= , 



(3a) 



(3b) 



(4) 



where D/rn? is the adatom diffusion coefficient in the 
velocity space in the absence of the StAW and 7 is the 



friction coefficient, given by 7 = ksT/D. Eqs. ((3)) re- 
lies on the high number of atoms in the substrate and 
its thermaUzation by an external macroscopic system in 
experimental conditionsJ^iH In this model, the correla- 
tion time Tc is of the order of the inverse of the Debye 
frequency of the substrate. The meaning of the friction 
coefficient appears clearly in the limit of vanishing corre- 
lation time : 



where 6 is the Dirac distribution, and thus 



,-,dx , , Te^o dx 



(5a) 



j(:,it),xit),t~t')-it')dt' ^=^^-it), (5b) 



dt 



dt 



the usual friction force proportional to the speed. In 
this limit, 771/7 is the relaxation time tr of the ad-atom 
dynamics in the absence of r.h.s. in Eq. ([T]). 

Eq. (H)) involves two different length-scales, the StAW 
wave-length A = 27r/fc and the periodicity a of the crys- 
talline potential; and three different time-scales, the cor- 
relation time Tc, the StAW period T = 2tt/uj and the 
relaxation time r^ of the ad-atom dynamics. In the fol- 
lowing analytical calculations, we assume that all these 
scales are separable and have the following properties: 



a< A 

Tc<^T «.TR 



(6a) 
(6b) 



Eq. (|6aP derives from the fact that we consider StAW 
wavelengths varying from few to hundreds of nanome- 
ters compared to the substrate lattice parameter of the 
order of 10^^ nm. Eq. (j6b[) is motivated by ours Molecu- 
lar Dynamic (MD) simulations of ad-atom diffusion on a 
substrate submitted to a nanometer wavelength StAW<^ 
The effective crystalline potential ($off) will be mod- 
elled by its fundamental Fourier component: 



$off(a;) = — --^cos(27ra;/a), 

ZTT 



(7) 



where (po is the amplitude of the corresponding force. 
Because of the difference of length scales (Eq. (|6ap ) , we 
assume for simplicity that the wavelength of the StAW 
force is an integer multiple of the lattice parameter A = 
77a. 



$off(a;) 



nk 



cos{nkx). 



(8) 



Within the presented model, Eq. ([T]) writes now: 



-\t-t'\ dx , ,, , , 
e^^—{t')dt' 



d^x 7 

df^ Tc J -00 dt' 

—ipo sm{nkx) + ^(t) + Fsaw sin(fca;) cos{ujt), (9) 

where ^ is the Gaussian noise defined by Eqs. ([3]). 

In a reduced dimensionless formulation Eq. ^ writes: 



X + — 

Tr 



' xdt 



with X — kx, i = ujt, 7 = ^/{muj), Tc — loTc, ^(t) ~ 
k^{t) / (muj'^) and F = kFsAw/i'>TT'^'^) ^^nd where x and 
X are the first and second derivatives of x relative to t 
respectively. 

Getting an explicit expression of the general solution 
of the non-linear stochastic intcgro-differential equation 
([9]) or (jlOp is unreachable at least for us and, anyways is 
not our goal. As stated before, what we are interested 
in, is to evidence the conditions for the appearance of 
a structuring effect on the ad-atom diffusion due to the 
effective StAW force and how the other forces can affect 
it. 

In Sect, ini we evidence the structuring effect of the 
StAW force, by considering Eq. PH)) in the long time 
(t 3> Tc) and length scales {x ^ a) limit without thermal 
fluctuations. In this limit, the periodic potential and the 
stochastic forces can be neglected and the retarded effects 
in the friction force too (limit Tc —s- 0). Eq. (fTO)) reduces 
then to: 



X + "fx — i^sin(i:) cos(t). 



(11) 



-ifo sin(77x) 4- ^(i) + Fsin(2:) cos(f), (10) 



Since Eq. (|lip is non-linear, we apply, in Sect. |TT] some 
standard tools of the non-linear physics to characterize its 
solutions. First, using the multiple time scale analysis^^ 
we will evidence the existence of an effective potential 
Ucs governing the diffusion in the long time scale limit. 
The general solutions will then be studied using the fixed 
points stability analysis. Then, it will be numerically in- 
tegrated and its solutions will be analyzed with the use 
of the Poincare sections of the phase diagram and the 
calculation of their Lyapounov exponents. The stochas- 
tic force S,{t) will then be reintroduced in Sect.|IlT]in the 
same Tc ^- limit. It will be shown that it mainly in- 
duces some fluctuations around the solutions of Eq. (|lip . 
In Sect. IIVI the retarded effects {tc 7^ 0) will be added, 
whereas the effects of the substrate effective cristalline 
potential will be reintroduced in Sect. |V]in the absence 
of retarded effects (limit Tc — > 0). In all the sections 
from mil to |V] , the structuring effect of the StAW force 
and its sensitivity to the other forces will be evidenced 
through the study of the position probability density of 
the ad-atom. Finally, in Sect. IVII the relevance of the 
analysis of Sect. |TT]on the solutions of the complete equa- 
tion (Eq. ([TU)) ) will be demonstrated. 



II. MAIN EFFECTS OF THE STANDING 
ACOUSTIC WAVE ON AD-ATOM DIFFUSION 



We study in this section the non-linear equation 
Eq. (jlip . Note that this equation is invariant under a 
simultaneous space and time translation: a; — ^ a; -I- tt and 
t — !■ < + TT. As a consequence, all the results concern- 
ing the structuring effects will be invariant under a space 
translation x ^- x + n. 



A. Multiple time scale analysis 

Eq. (jlip is a non-linear deterministic equation that im- 
plies two different typical time scales: the dimensionless 
relaxation time f^ = ujtji (depending on 7) of the ad- 
atom dynamics and the period of the StAW force. To 
manage these time scales, this equation is first analyzed 
using the multiple scale method.— This method, due to 
the extra degrees of freedom it introduces, allows to re- 
move the secular divergencies that can arise in a standard 
perturbation approach. Note that Eq. (fTT]) without fric- 
tion has already been studied^^ in a different framework. 

We consider the limit e = l/w -^ (keeping constant 
Cs = u)/k, the sound speed in the substrate). Eq. (fTTj) 
reads then: 



X + ejx = eFsin(a;) cos(t), 



(12) 



where 7 = 7/m and F_ = FsAw/i'mCs) are order 
quantities (0(e°)) . We seek an approximate solution 
of Eq. p^ using the following expansion: 



x(i,e) = xo{io,ii) + exi(io,ii) 



(13) 



involving the two time scales: to ^ t and <i = ei. Sub- 
stituting Eq. (fT3|) into Eq. (|12p and identifying terms of 
the same order in e, we get the following equations: 



DgXo = 0, 



(14a) 

DIxi — —2DoDiXo - jP^xq -I- Fcos(to) sin(2:o), 

(14b) 

Dli2 = -2DoDi5:i - DJia - i{Dq3:i + Diio) 

+F_xi cos(io) cos(xo), (14c) 

where the operator D„ designs d/din with n e {0, 1}. 
The solution of Eq. (|14ap reads: 



xo=A{h)io + B{h). 



(15) 



The first term of Eq. p5|) . a secular term that diverges 
with ip: is removed by setting A{ii) — 0, leading to: 



io = £o(^~i)- 



(16) 



Substituting Eq. ([TC]) into Eq. (J14b|) . the particular 
solution x\ writes: 



ii(ioJi) == -£cos(to)sin(xo). 



(17) 



Using these expressions of xq (Eq. ((T5)) ) and x\ (Eq. ([T7)) ) 
in Eq. p4cp . we write the solubility condition of this 
equation by the elimination of the secular term: 



F' 



D'(xo + 7D1X0 = -^ cos(xo) sin(£o), (18) 

a differential equation governing the solution xq on the 
time scale ii. To give a physical meaning to Eq. ([T51) . 



we note that using Eqs. ([16]) and (H?]) in Eq. (HS]), the 
solution of Eq. (|lip writes to O(e^) order: 



x{i) = .To(ti) - Fcos{io) sin(£o(ii))- 



(19) 



Therefore, X{ii) the average value of x{i) over a StAW 
period, writes to O(e^) order: 



X{h) = (i(t)} = — 

ZTT 



to+27r 



x{i)difv£o{ii). (20) 



Substituting xo(ii) by X{ii) in Eq. ([T8)) an going back 
to the t variable, we obtain the long time evolution equa- 
tion of X{i) ~ {S:{i)), the mean value of x over a StAW 
period: 



d^X{t) ^dX{i)) 



df^ 



dt 



■sin(2X(t)). (21) 



Note that the scheme provided by Landau and Lifshiftii 
that develops the variable x{t) in Eq. ([TT|) as a sum of a 
slowly varying function X(t) and a quickly varying func- 
tion C(t) yields a similar result 1^ 

The multiple scale method has allowed us to trans- 
form the non-autonomous Eq. (jlip into an autonomous 
equation (Eq. (|2ip ) on a longer time scale. This equation 
describes the motion of the ad-atom on a mesoscopic time 
scale, long compared to the period of the StAW, but small 
compared to the relaxation time fn of the ad-atom dy- 
namics. The StAW force acting at the mesoscopic time- 
scale derives from the effective potential UcS'- 



U,s{X) = ^W{X). 



(22) 



This potential, is periodic with period tt and minima at 
X = 0[7r]. Whatever the initial conditions the long time 
evolution described by Eq. (PT|) will be a dampted evo- 
lution towards one of the minima of C/eff . Within the de- 
scribed approximation (Eq. ((2T|) ). the StAW leads then 
to a self-organization of the ad-atoms diffusion into a pe- 
riodic array with period half the period of the StAW. In 
the following, UeS will appear to be an essential tool to 
interpret the ad-atom trajectories and the structuring ef- 
fect of the StAW. Note that the (stable) fixed points of 
Eq. (|2ip which are the minima of C/eg, are also those of 
Eq. dn]). 

These results are approximated results, we must now 
come back to Eq. pT|) to test their relevance in the gen- 
eral case. We will start with the study of the stability 
of the fixed points {{x,x) ~ (0[7r],0)) of Eq. ([TT]) which, 
as we will see now, can differ from that of Eq. ([21]) for 
certain values of the parameters of the equation. 



B. Fixed points stability 

As already mentionned, Eq. (jlip is invariant under 
simple spatial and time translations so that all the fixed 



(?- 




FIG. 1. Stability (Blank) (SD) and instability (dashed) (ID) 
domains of the fixed point (0,0) of Eq. ((TTJ in the (7,-?") pa- 
rameter plane (PP). Red crosses design the couples of param- 
eters chosen in the numerical studies of Sects. HIl IIIII and IIVI 



points {x,x) = (0[7r],0) are equivalent. We hence reduce 
our stability analysis to one fixed point: {x,x) = (0,0). 
The linerarized version of Eq. (|lip around this fixed point 
reads: 



V = A{i)V{i) 



(23) 



with V = 



and A{i) 



1 

Fcos{i) -7 



We define Rr , the propagator of Eq. (|23] 



to 



RL 



V{to) ^ Vit) 



(24) 



Since A{i) = A{i + 2tt) in Eq. 1^, the Floquet theory 
provides the stability of the fixed point {x,x) = (0,0) 
of Eq. (|TT|) from the eigenvalues of Rq^: the fixed point 
is stable if all the eigenvalues of Rq"^ arc inside the unit 
circle of the complex plane. 

The propagator R^'^ is calculated by numerically in- 
tegrating Eq. (|23p using a fourth order Rungc-Kutta 
method and is then diagonalized. Fig. [T] reports the sta- 
bility diagram of the fixed point (x, x) = (0, 0) in the 
(7, F) parameter plane (PP). Fig. [T] reveals the existence 
of unstable domains which questions the validity of the 
self-organization effect evidenced for any (7,i^) parame- 
ters in the preceding section. This is not contradictory 
since the Floquet analysis is exact concerning the fixed 
points stability, while the approximated multiple scale 
analysis relies on the smallness of the amplitude of the 
fast variations (StAW period time scale) compared to the 
slow variations {fn time scale) of x. A condition which 



will be fulfilled as far as the strength of the friction force 
(7) remains large enough compared to the strength of the 
effective StAW force (F), which explains the global sepa- 
ration between the stability (upper left triangle) domain 
(SD) and the instability (lower right triangle) domain 
(ID) in the PP. In addition, in absence of friction 7 = 0, 
Eq. (|23p is equivalent to a Mathieu equation for a para- 
metric oscillator with a null eigen frequency. The stabil- 
ity diagram of the Mathieu equation can be for instance 
found in Fig. 2 of Ref . [l8t it evidences some unstable solu- 
tion regions, corresponding to the parametric resonances, 
separated by very tiny stable regions for an oscillator 
with a null eigen frequency. The large separated unsta- 
ble subdomains (tongue) in Fig. [T] hence correspond to 
the parametric resonances that disappear when increas- 
ing the dissipation 7. 



C. Numerical study 

In order to check the self-organisation ability of the 
effective StAW force in the ID of the PP (Fig. [T]), we 
need to turn to a precise examination of the trajectories. 

An analytical resolution of Eq. ([TT|) is out of scope, 
so we will solve it numerically using the forth order 
Runge Kutta method for different values of the parame- 
ters (7,i^). We will also look at the trajectories for pa- 
rameters in the SD to give a complete scope of the differ- 
ent exact behaviors. In order to enlighten some specifici- 
ties of the trajectories, the resolution will be performed 
for different initial conditions {x{0),x{0)) (IC). From our 
analysis of the stability diagram in the preceding sec- 
tion, we expect to observe specific behaviors associated 
to increasing values of F at constant 7, going from one 
stability domain to the next one through the midway in- 
stability domain. This is what we observed in our rather 
extended exploration of the PP and IC spaces. From this 
exploration, we identified four categories of trajectories: 
one corresponding to the SD; the three others to the ID, 
one in its core and the two others in the vicinity of its 
frontiers with its two neighboring SDs. We did not ex- 
amine the very peculiar case 7 = 0, which corresponds 
to undamped trajectories, since it is not relevant for the 
ad-atom diffusion on a substrate, and since it has been 
studied previouslyJ^ Of course, we cannot absolutely ex- 
clude the possibility to have missed some specific behav- 
iors, even if we consider it as highly improbable. In all 
the rest of the manuscript, we will focus on the solutions 
for the 2.0 constant 7 value, since it provides a represen- 
tative sample of the behaviors we have exhibited. The 
results are presented for increasing values ofF (1.5, 3.15, 
4.5 and 4.96) through the first instability subdomain (red 
crosses in Fig. [1]). In each case, the trajectory is related 
to the effective potential Ues{X) ( blue solid line in Fig. [5] 
and Eq. (j22p ) in order to evidence any self-organization 
behavior. 



1. Converging Trajectories 

The trajectories have typicahy an osciUating behavior 
whose amphtude and off-set are decreasing. They can 
be somewhat different at the very beginning, depending 
on the sign of the 5:(0)i(0) product but they have the 
same character at a longer timescale. Fig. [2^ reports the 
solution for initial conditions (1.0,0.0) corresponding to 
an initial position in the Uefi{X) potential valley associ- 
ated to its 0.0 minimum. We choose an initial position 
(1.0) rather away from the minimum of Ucii{X) to best 
evidence the amplitude and off-set decreases. The initial 
speed has been fixed to zero in order to only present the 
characteristic time behavior. 

The trajectories converge towards the fixed point 
(0, 0). Depending on the initial conditions, the trajectory 
can eventually escape from the 0.0 minimum to an ad- 
jacent minima. In that case the trajectory will converge 
to the corresponding Ues{X) minima (fixed point). The 
several time scales mentioned in Sect. Ill Al are clearly vis- 
ible: fast oscillations at a 2tt period (StAW period) whose 
amplitude slowly decays on the tj^ time scale. The pe- 
riod of the oscillations can be one of the harmonics of the 
StAW period. For example, in the vicinity of the stabil- 
ity/instability frontier (7, -F) = (2.0,3.028), we observe 
the second harmonic period (47r). It appears then that, 
in the SD, the solutions of Eq. ([TT|) can be adequately 
estimated using a multiple time scale analysis, the long 
time behavior being well described by Eq. (j2ip . 



Effective polcnlial (iirb. unil) 
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FIG. 2. Black and red curves: x{t) solutions of Eq. (|lip 
for 7 = 2.0. Stability domain (SD): a) converging trajec- 
tories [F — 1.5). Instability domain (ID): b) periodic tra- 
jectories {F — 3.15), c) chaotic trajectories {F = 4.5) and 
d) unbounded trajectories {F — 4.96). Initial conditions are 
reported in the legend. Blue curves: UcsiX), the effective 
potential (abscissa) in arbitrary units versus position X (or- 
dinate). 



2. Periodic trajectories 

In the ID, in the vicinity of the first frontier be- 
tween the SD and the ID, the unstable character of the 
fixed point results in an amplification of the oscillations 
(Fig.[2lD). This amplification occurs up to an upper limit 
fixed by the non-linearities of the sinus in Eq. ((TT|) . At 
this time the trajectory becomes periodic, revealing the 
existence of a limit cycle, the attractor of the system 
in the phase space. The IC of the presented trajectory 
(i(0),a;(0)) = (0.1,0.0) correspond to a position very 
close to the unstable fixed point to evidence the amplifi- 
cation of the fast oscillations (transient regime) leading 
to the periodic trajectory. The period of oscillation is 47r, 
i.e. twice the excitation period. Like in the converging 
trajectories, this period depends on the peculiar choice 
of the parameters (see Sect. HID)) . 

The trajectory of Fig. [JId oscillates in its initial effec- 
tive potential basin. However, increasing F, results in an 
increase of the upper limit. This limit may overcome the 
maxima of Ucff{X) so that the trajectory may extend on 
the two neighboring effective potential basins: however, 
the average position of the adatom always belong to the 
same potential basin. Ucs{X), if relevant to characterize 
the trajectory of the adatom, does not capture all the 
information contained in the phase space diagram. With 
this kind of trajectories the self-organisation character 



of the effective StAW force is not apparent but will be 
demonstrated in Sect. Illll studving the distribution of the 
ad-atom x-position. 



3. Chaotic trajectories 

In order to evidence the chaotic character of some tra- 
jectories in the core of the ID region {F = 4.5) we re- 
port on Fig. [2t the solutions of Eq. PT|) for two sets of 
very closed I.C.s (0.10,0.0) and (0.12,0.0). Both trajec- 
tories visit different effective potential basins following 
an apparently erratic motion. The two solutions, diverge 
quickly from each other, despite their very close ICs, sug- 
gesting the chaotic character of these trajectories. To ex- 
amine this chaotic behavior, we calculate the Lyapounov 
coefficients (A.;) of the autonomous system associated to 
Eq. (fTTj) . For the first trajectory, using the algorithm 
provided by Wolf ct ali^^ the highest Lyapounov coef- 
ficient Al is found to be positive (0.278) (binary base), 
an unambiguous evidence of the chaotic character of this 
solution. 

The transition to the chaotic behavior will be discussed 
in a separate section fSect. [irD|) . At first sight, it seems 
that the self-organization character is lost, but as we will 
see later (Sect. HIT)) studying the position probability den- 
sity of the ad-atom, it actually is preserved. 



4- Unbounded trajectories 

We look now, in the ID domain, at a charaeteristic 
trajectory in the vicinity of the second frontier between 
the ID and the SD (F = 4.96). The IC of the trajectory 
reported in Fig. [5Ji are the same as those in the two pre- 
ceding cases, (0.1,0.0). After a transient period where 
the particle stays in its original potential well, it leaves 
it, without being captured by any other potential well: 
the trajectory is unbounded. As one can see, from this 
point the trajectory is roughly linear and thus does not 
present any visible chaotic character. This is confirmed 
by the calculation of its Lyapounov exponents that are 
all negative (or null). After the transient period, the 
trajectory is monotonously increasing with a staircase 
character. Changing the IC can lead to a monotonously 
decreasing trajectory with a symmetric staircase charac- 
ter. These trajectories are analogous to the rotations of 
a pendulum about its pivot in the clockwise or counter- 
clockwise directions and have been already evidenced in 
the phase space (a;(0),x(0)) in the absence of dissipation 
7 = 0J^ 

The final remark of the preceding section still applies 
here: contrary to the appearances, the self-organization 
character is preserved. It infers from the staircase char- 
acter of the trajectory which corresponds to longer res- 
idence times in the potential valleys than the transi- 
tion times between valleys. The self-organization char- 
acter will be evidenced for all the types of trajectories in 
Sect lllli where the probability distribution will be stud- 
ied in the presence of the fluctuating force (vanishing 
fluctuating force case). 



D. Transition to Chaos 

We will now characterize the domains of existence, in 
the PP plane, of the four observed types of trajectories, 
through the study of the Poincare section {x = 0) of 
the solutions of Eq. ((TT|) at times larger than the tran- 
sient initial period. As in the preceding section, wc ex- 
plore the PP at the constant 7 = 2 value. We increase 
F starting with a value F — 1.5 in the flrst SD up to 
F = 6.84 in the heart of the second ID, going through 
the flrst ID in between (Fig. [T]). The critical values of F 
separating different behaviors have been calculated using 
the shooting and continuation methods2S and are given 
here with a precision of 0.001. Note that these F criti- 
cal values depend on the 7 value. The calculations have 
been performed for a wide range of ICs. Whatever the 
parameters, we found at most two types of asymptotic 
trajectories presented in Fig. [5^ (black and red dots) de- 
pending on the ICs. The presented points correspond 
to (0.1,0.0) (black dots) and (-0.1,0.0) (red dots) ICs. 
Al, the greatest non-null Lyapounov coefficient (binary 
base) of the autonomous system associated to Eq. ([TT|) , 
has been calculated also to characterize the chaotic or not 
character of the trajectories^^ Since the results are in- 
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FIG. 3. a) Poincare section of the phase space omitting tran- 
sient regime as a function of F with 7 — 2.0 from trajecto- 
ries calculated by numerical resolution of Eq. (|ll|l with initial 
conditions (0.1,0.0) (black) and (-0.1,0.0) (red). Blue curve: 
Ucs{X), the effective potential (abscissa) in arbitrary units 
versus position X (ordinate) b) lambdaL the greatest non-null 
Lyapounov exponent (binary base) of Eq. (|lip as a function 
of F for initial conditions (0.1,0.0). The solid horizontal line 
[lambdaL = 0) is a guide to the eyes. 



dependent of the ICs, we present in Fig. ^p its evolution 
with F for the (0.1,0.0) IC only. 

From Fig. [3] we see that in the first SD i.e. F < Fi, 
with Fi = 3.028 the first SD/ID limit, the asymptotic 
solution of Eq. (fTTj) is the stable fixed point (0, 0) what- 
ever the ICs (sec. lIlCTI) . This result is coherent with the 
values of Xl < 0. 

At Fi, Xl goes to for the first time and we enter the 
first ID, the system undergoes a Hopf bifurcation towards 
an unique limit cycle (with twice the St AW period), the 
attractor of the flow (identical black and red points) in 
agreement with the results of sec. Ill C 21 More precisely, 
due to the translation invariancc of Eq. (|lip . there is 
actually an inflnite number of limit cycles, attractors of 
the system, one per flxed point (0[7r], 0). This domain of 
periodic asymptotic trajectories extends from Fi to F2 — 
3.654. Beyond, for F > F2, Xl becomes positive and 
thus the trajectories chaotic. The flrst periodic-chaotic 
transition is thus at F = F2. 

In this periodic domaini^i < F < F2, a.t F — 3.552, 
Xl goes to and the limit cycle splits into two limit cy- 
cles (separation of the black and red points) without any 
period change. The system has then two attractors per 
flxed point. This doubling of the number of attractors, 
precedes for each attractor a cascade of period doublings 
at increasingly close values of F, each of them being as- 
sociated to the cancellation of Xl ■ We have been able to 



observe 4 period doublings (at 3.632,3.650,3.653,3.654), 
though we did not try to optimize this number. Each 
period doubhng corresponds to the apparition of the cor- 
responding peak in the Fourier spectrum. The transition 
to chaos proceeds when the spectrum becomes continu- 
ous at F2. 

Between F2 and F3 = 4.890, there is an alternation 
between large domains of chaotic trajectories and many 
very small domains of periodic or unbounded asymptotic 
solutions (with Al < 0) . The chaotic trajectories are 
characterized by positive values of Xl and Poincare sec- 
tions containing an infinite number of points in the limit 
of an infinite trajectory. In the very small domains where 
Xl becomes negative, the solutions arc either asymp- 
totically periodic with, in the Poincare section, a re- 
duced number of points, or unbounded with an absence 
of points in the x-interval chosen for Fig. [3^. 

At F3 — 4.890, Xl goes again to and there is a transi- 
tion towards unbounded asymptotic trajectories {Xl < 
and absence of points in the x-interval of the Poincare 
section Fig. [3K)- This unbounded domain extends up to 
the end of the first ID at F = i^4 = 5.090 where again 

For F4 < F < F^, where F^ = 6.183 corresponds to the 
second SD/ID limit, the asymptotic solution of Eq. ([TT|) is 
again the stable fixed point (0, 0) with Xl < 0, whatever 
thelCs fsec HTCTI) . 

At F5 we enter the second ID and wc observe a be- 
havior very similar to the one in the first ID. The main 
difference is that the first periodic domain starts directly 
with two limit cycles, then there is also a cascade of pe- 
riod doublings, the first one at F = 6.522, leading also 
to a region of alternation of wide chaotic domains and 
small periodic or unbounded domains. 

The system described by Eq. (|lip is very rich and com- 
plex. The apparition of chaotic solutions is actually not 
surprising: imposing a standing acoustic wave involves 
the interaction of two harmonics in a parametric-like ex- 
citation. The equation (|11|) without dissipation has been 
studied in order to investigate regimes where resonances 
induced by both harmonics interact J^ The kicked rotor- 
is also very similar to our system excepted that an infin- 
ity of harmonics are involved in the excitation, so that 
their resonances can interact leading to chaotic solutions. 
Finally let's mention the work of Van Dooren^S who stud- 
ied the dynamics of a pendulum with a forced sinusoidal 
horizontal support motion: in the absence of gravity, this 
system reduces to our undamped system (7 = 0). 

From this numerical study, it appears that the trajec- 
tory of an ad-atom on a substrate submitted to a StAW 
can be of four different types, depending on the parame- 
ters (7, F). While increasing values of F at constant 7, 
the domains corresponding to the different kinds of tra- 
jectories are successively i/ converging trajectories in the 
SD domains, followed in the ID domains by ii/ a periodic 
domain, then iii/ an alternation of mainly chaotic solu- 
tions together with small periodic or unbounded domains 



and ends up with iv/ an unbounded domain, and so on 
when entering the next SD. While the structuring effect 
of the StAW is obvious in the case of converging trajec- 
tories and to a less extend periodic trajectories, since it 
drives the particle into given regions in the configuration 
space, it is less obvious for the other cases, in particular 
for the chaotic domains. 



III. STOCHASTIC FORCES 

In this section, we reintroduce the Gaussian stochastic 
force ^(i) in Eq. pT|) in the Tc — ;• hmit: 

x{i) + ^x{t) = ^{t) + Fcos(i) sin(i), (25) 

with 



l(i))=0. 



(26a) 



(26b) 



where D ^ D/{ujm'^cl) (Eqs. © and ©). Note that 
D and 7 are not independent due to the fluctuation- 
dissipation theorem (Eq ^): 
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(27) 



We will use in this section the same 2.0 fixed value of 7 
and the same four values of F as in the previous sec- 
tion, corresponding to the four exhibited types of so- 
lutions. Three values of D (0.0001, 0.01 and 0.1), or 
equivalently of temperatures for given substrate (cg) and 
ad-atom (m), will be investigated, corresponding to the 
quasi-absence of the stochastic force (Sect. HIl Eq. (fTTjl ). 
a medium and a strong stochastic force respectively. 

Due to the stochastic character of Eq. (|25p we adopt 
here a statistical point of view. It appears, as demon- 
strated in this section, that on a longtime scale com- 
pared to the relaxation time 77^, the memory of the IC 
is lost and the system is in a quasi-steady state: P{x), 
the distribution of the ad-atom x-position is mainly time- 
independent (rigorously, it involves a tiny periodic con- 
tribution at the StAW frequency) . On a timescale larger 
than the StAW period, the structuring effect of the StAW 
will be revealed through the correlation between the x 
variations of P{x) and those of the periodic effective po- 
tential [/off- Physically, P{x) will point out the preferen- 
tial sites where the ad-atom spends most of its time. 

Assuming the ergodicity of the system, P{x) can be ob- 
tained from a single long trajectory of one particle (after 
elimination of the initial transient period, with whatever 
IC) or from a set of trajectories. In addition, due to the 
translational invariance x^-x + TT,t^t + TToi Eq. pip . 
P{x) is expected to be 7r-periodic. Wc calculate P{x) 
from a number of trajectories obtained for different ICs 



and realizations of £,{t): 100 trajectories of t = 10000 
time units each, The ICs are taken at random in the 
a; = effective potential valley with x ~ 0. The points 
of the trajectory outside the initial valley are translated 
back to the i = valley coherently with the translational 
invariance of Eq. (fTT]) . P{x) is then estimated from the 
histogram of the ad-atom position of these trajectories 
and the following normalisation condition: 



7r/2 



P{x)dx = 1. 



(28) 



-7r/2 



The results are reported in Fig. |4] over half the StAW 
wavelength, i.e. a period of UcS together with the 
effective potential UcB- 



1. Converging trajectories 

F = 1.5, Fig. Wn. At a very low diffusion coefficient 
value [D = 0.0001), P{x) is strongly peaked at a; = 0: in 
the quasi-absencc of fluctuating force, the trajectory still 
converges to the minimum of the effective potential. The 
tiny stochastic force induces small fluctuations of the po- 
sition in the vicinity of the minimum. This fluctuations 
are not strong enough to induce a transition to an ad- 
jacent valley (inter- valley transition) on the simulation 
time scale (P(±7r/2) = at the maxima of C/eff)- At 
D = 0.01, the width of the peak centered on the fixed 
point of Eq. pT|) has not sufficiently increased to induce 
a significant intcr-valley transition, whereas it does at 
i) = 0.1 (P(±7r/2) ~ 15%): the ad-atom diffuses from a 
basin of attraction to a neighboring one. The converg- 
ing character of the trajectory is lost. However the self- 
organization is preserved through the peaked character 
of P{x) centered on the minima of UcS- 



2. Periodic trajectories 

F ~ 3.15 - Fig. ^p. Here also at very low diffusion co- 
efficient D = 0.0001, the periodic character of the trajec- 
tory is roughly preserved with small fluctuations around 
the initial periodic trajectory in the absence of ^; and the 
particle visits a wide region of a basin of attraction of the 
effective potential. The peaks and features observed on 
the plot of P{x) for D — 0.0001 arc due to the speciflc 
shape of the limit cycle (or trajectory shown in Fig. [^b). 
Increasing the diffusion coefficient D induces some fiuc- 
tuations around this limit cycle. They can even activate 
the crossing of the effective potential barriers (on the 
simulation time scale), clearly evidenced for D = 0.1 by 
the significant value of P{ti/2). Increasing the diffusion 
coefficient smooths the structural role of the StAW: the 
stochastic fluctuations give rise to a wide Gaussian-like 
distribution centered on the effective potential minimum 
and whose width increases with the diffusion coefficient. 



3. Chaotic trajectories 

F ~ 4.5 - Fig. [^. In the quasi-absencc of stochastic 
force [D — 0.0001) the chaotic character of the trajec- 
tory leads to a P{x) distribution correlated to UcS'- even 
if the ad-atom is not trapped in a given potential val- 
ley (P(7r/2) 7^ and Fig. [2t), P{x) has a pronounced 
maxima at the minima of C/cff. Increasing the diffusion 
coefficient yields the same qualitative observations as in 
the previous periodic trajectory case. 



4. Unbounded trajectories 

F = 4.96 - Fig. g^. At the very low D = 0.0001 value, 
from the staircase character of the trajectory presented 
in Fig.[2}l one expects a non uniform P{x) with a marked 
peak inside each potential well. A precise examination of 
the trajectory in Fig. [2Ji reveals that the plateaux are at 
X positions slightly larger than the center of the potential 
wells. Such a trajectory contributes then to the ii > 
peak observed on Fig. |3Ji (top). The second peak at the 
symmetric —xi position results from mean trajectories 
towards a; < positions. The ii > non centered posi- 
tion of the distribution associated to trajectories towards 
increasing mean x values results from the definite direc- 
tion a; > or X < of the observed trajectories (a; > in 
Fig. [2}1) . Increasing the diffusion coefficient results in an 
increased width of each peak, leading to a unique central 
peak for D = 0.1. As in the three preceding cases, the 
structuring effect of the StAW is also evidenced in that 
case whatever the strength of the stochastic force in the 
studied range. 

As a conclusion, we would like to emphasize that for 
all the trajectories types, and as soon as the stochastic 
force is significant, the distribution P{x) is a Gaussian- 
like distribution centered on the effective potential min- 
imum. Therefore, the StAW has a structuring effect on 
the diffusion of the particle: it induces preferential sites 
in which the particle spends more time. These sites pre- 
cisely correspond to the minima of the effective poten- 
tial evidenced in Sect. [Ill The stochastic force essentially 
counterbalances the structurating role of the StAW by 
smoothing the distribution P{x). 



IV. RETARDED EFFECTS 

We have neglected the correlation time Tc in Sects. [ll] 
and mil In surface diffusion problems, such an approxi- 
mation is valid for heavy ad-atomsi^ which is a very pecu- 
liar case. In this section, we investigate the effect of a sig- 
nificant value of the correlation time Tc compared to the 
StAW period. The correlation time Tc is involved in the 
friction term through the memory kernel (Eq. ^) and in 
the auto-correlation function of the stochastic force ^(f) 
(Eq. ()3bp ). Reintroducing the retarded effects Eqs. (|25p . 
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FIG. 4. Black, red and green curves : Histograms of P{x) 
for different values of the diffusion coefficient D: 0.0001, 0.01 
and 0.1 respectively, with 7 — 2.0 and different F values, 
i.e. types of trajectories as defined in sec. Ill Cl (a) converging 
(F = 1.5), (b) periodic {F = 3.15), (c) chaotic [F = 4.5) and 
(d) unbounded (_F = 4.96). Blue curve: Uch{X), the effective 
potential in arbitrary units. 



P6ap and (j26bp write now: 



i{i) + ^ I .i(<')e-(*-*'^/"^dt' = F cos(t) sin(.i) + |(t), 



l(i))=0, 
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(29) 
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with the same Eq. ^7} between D and 7. 

Let's first estimate a physical range for the time Tc = 
ujTc- On the one hand, in the case of the diffusion of 
an ad-atom, 1/tc is roughly of the order of the De- 
bye frequency f]j, corresponding to the maximum fre- 
quency of atomic vibrations in the crystalline substrate, 
i.e. fu « 10^'^ Hz for common crystals, and consequently 
Tc « 0.1 ps. On the other hand, as mentioned in the in- 
troduction, the StAW wavelengths of interest vary from 
few to hundreds of nanometers, i.e. typically from 5 nm 
to 1 ^m. With typical sound speeds in solids around 3000 
m.s^^. the StAW frequency uj lies in the range 1.25 10^° 
- 2.5 10^^ rad.s"^. Consequently, fc will be in the range 
0.002 - 0.40. 
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FIG. 5. Converging trajectories x{t) solutions of Eq. (|29|) . 
(7, F) — (2.0,1.5), with a small diffusion coefficient {D = 
0.0001) and initial conditions {x,x) = (1.0,0.0), for three dif- 
ferent values of fc: 0.001 (black), 0.3 (red), 0.4 (green). Blue 
curve : Ucff{X), the effective potential (abscissa) in arbitrary 
units versus position X (ordinate). 



Three typical values of fc, 0.001, 0.3 and 0.4 will 
be used to investigate the effect of Tc in the two ex- 
treme cases: quasi absence {D = 0.0001) and significant 
(D = 0.01) fluctuations, two values of D already used 
in the preceding sections. Due to the separation of the 
different time scales, we expect the solutions of Eq. ([^^ 
with a low fc value (0.001) to be very similar to solutions 
of Eq. ((25|) . Only values of fc non negligible compared 
to the StAW period are expected to produce solutions of 
Eq. ([^ significantly different from those of Eq. (^5]) . We 
reduce our study to the converging trajectory case pre- 
viously studied (7 = 2, F = 1.5), with the same (1.0,0.0) 
fC as in Sect. IIIC 11 since this type of trajectory is the 
most favorable for self-organisation and thus will be usu- 
ally preferred in any application. 

Eq. (p9|) is numerically solved using a Leap-frog 
algorithm)^ ^(i) values satisfying Eqs. (I30al) and ()30bp 
arc generated with the algorithm of Ref. |23| while the in- 
tegral of the friction term including the memory kernel 
are calculated using the algorithm given in Ref. Il3l 

The solutions for D — 0.0001 are reported on Fig. [S] 
As expected, when the fluctuations arc negligible, what- 
ever the correlation time, the solutions of Eq. ([TT|) are 
qualitatively unchanged and thus Tc has a negligible in- 
fluence on the structuring effect. The solutions are still 
oscillating functions at the StAW time scale and they 
still converge to the fixed point. The main effect of in- 
creasing values of Tc is an increase of the amplitude of 
the oscillations at all timescales, and correlatively of the 
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FIG. 6. (a) : converging trajectories x{t) solutions of Eq. (|29|l . 
{■y,F) = (2.0, 1.5), with a diffusion coefficient D = 0.01 and 
initial conditions {x,x) = (1.0,0.0), for three different val- 
ues of fc'. 0.001 (black), 0.3 (red), 0.4 (green). Blue curve : 
Uch{X), the effective potential (abscissa) in arbitrary units 
versus position X (ordinate), (b): Corresponding histograms 
of P{x) from 100 trajectories of 10 000 time units with initial 
conditions evenly distributed between — tt and tt (same color 
code as in (a)). 



relaxation time tji of the ad-atoni dynamics. Such be- 
havior has already been mentioned and explained in the 
literature: the velocity autocorrelation function for the 
ad-atom shows both a ballistic and a diffusive regime, 
the width, roughly ffl of the balistic regime, is related to 

At odds, when the stochastic force is large (D = 0.01), 
the possibility for the ad-atom to cross an effective poten- 



tial barrier on the simulation time scale increases with Tc'. 
we report in Fig. [S^ the trajectories of the ad-atom as a 
function of time for the three different mentioned values 
of fc- Increasing fc, and hence the ballistic regime, gives 
more importance to the very unlikely extreme values of 
^{t) and thus results in a higher effective diffusion coeffi- 
cient (not shown). Fig.|6lD reports the P{x) distributions: 
fc slightly affects the histogram P in the vicinity of the 
minimum of [/off- The effect is more pronounced on its 
wings: they increase significantly leading to non zero P 
values at the Ucs maxima for fc — 0.3 and 0.4, coher- 
ently with the possibility for the ad-atom to escape from 
its original potential valley evidenced in Fig. |B^). Nev- 
ertheless, no matter the value of fc, the shape of P{x) is 
still gaussian-like evidencing the structuring effect of the 
StAW 

We can thus conclude that if the retarded effects quan- 
titatively modify the trajectories, they weakly affect the 
structuring effect induced by the StAW. 



V. EFFECTIVE CRYSTALLINE POTENTIAL 

In this section, we consider the additional effect of the 
effective crystalline potential $off(i) on the motion of 
the ad-atom, in the same negligible correlation time limit 
(tc -^ 0). Eq. (Uni) writes then: 



x{t) + 7x(i) — ~(/?o sin{nx) + F cos{t) sin(.f) 



(31) 



The calculations will be performed with n = 24, a value 
comparable to the ones we used in our MD simulationsi^ 

First, we study the modifications induced by ^off on 
the fixed-point stability diagram described in Sect. |TT] 
- Fig. [T] in absence of stochastic forces. Note that in 
Eq. (|3T|) both the crystalline and StAW forces cancel for 
X = 0. We could easily imagine that a dephasing of $off 
compared to the StAW will shift the fixed points or yield 
to the absence of fixed points. However, since the StAW 
wavelength is large compared to the lattice parameter, 
we do not expect such a dephasing to qualitatively mod- 
ify the ad-atom trajectories, especially in the presence 
of the stochastic force. Fig. [7] reports, as in Fig. [I] the 
stability diagram in the {'j^F) PP for increasing values 
of ipo = 0.0, 0.05, 0.15, 0.25 and 0.35. The main effect o_f 
$off is to shift the instability domains towards higher F 
values at constant 7. The effective crystalline potential 
has then a stabilizing effect on the trajectories. 

To investigate the structuring effect of the StAW, as 
in the preceding section, we will study the effect of $cff 
in the presence of the stochastic force only for one couple 
(7 = 2,i^ = 1.5), corresponding to converging trajec- 
tories in the absence of crystalline potential. Here we 
are mainly interested on what happens when <I>off is sub- 
stantially higher than the effective StAW potential Ucff 
and is thus potentially able to challenge or to overcome 
the structuring effect of the StAW: actually, the StAW 
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FIG. 7. Stability (blank) and instability (dashed) domains of 
the fixed point (0,0) of Eq. (|31|l in the (7,-^) parameter plane 
(PP) for increasing values of the effective cristalline force (/?o • 
From top to bottom : ipo = 0.00 (same as Fig. [TJ, 0.05, 0.15, 
0.25 0.35. 



force becomes then a second order effect compared to 
the crystalline potential. In order to induce an efficient 
diffusion across the effective crystalline potential barrier 
on the simulation time scale, we will use the strongest 
stochastic force {D ~ 0.1) previously used. 

Solving numerically Eq. ([3T|). we calculate the his- 
togram P{x) from 100 trajectories of 50000 time units 
each with different ICs. Fig. [5] reports the distribution 
P{x) (normalized following Eq. (pS)) ) as a function of x 
for increasing values of 1^0 = 0.7, 1.5, 3.0 and 5.0, starting 
at a lower value than F up to 3.3 times F. Comparing to 
Fig El we see that the Gaussian like curves present now a 
structuration at the $off length scale. There are now two 
length scales: 1) a short wavelength oscillation (lenght 
scale a) due to 'I'off with the local maxima of P{x) at the 
minima of $eff ■ 2) a slow variation of the amplitude of 
the local maxima corresponding to the previous curves 
of Fig |6] and thus to the St AW with the maxima of P{x) 
at the minima of C/off- Unexpectedly, from Fig. |S1 the 
structuring effect seems to increase with the increasing 
of the effective crystalline potential. However, this effect 
is essentially due to the normalization: fitting the set of 
maxima of P{x) by a Gaussian curve leads to approxi- 
mately the same Gaussian width for all the values of (pQ. 

Hence, the structuring effect is weakly affect by the val- 
ues (even for significant values) of ipo ■ Note the very good 
qualitative agreement between the distribution P{x) in 
Fig.[5]and the results reported in our Molecular dynamics 
simulations^. 
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FIG. 8. Histogram of P{x) calculated from converging trajec- 
tories ((7,F) = (2.0, 1.5)) solutions of Eq. (|3TJ (100 trajecto- 
ries of 200000 time units each), with _D = 0.1 and increasing 
values of ipo from top to bottom: ifo = 0.7, 1.5, 3, 5. 



VI. EFFECTIVE POTENTIAL AND 
PROBABILITY DISTRIBUTION 

To quantify the structuring effect of the St AW, we note 
that the structuration can be characterized by an effec- 
tive energy difference AE^^yj at the mesoscopic scale be- 
tween the minima and maxima of the effective potential 
deduced from the P{x) curves of Fig. [U 






fcsTln [max(Pmax)/niin(Pmax)] (32) 



with Pmax the ensemble of local maxima of P(i;) and max 
(min) the ensemble maximum(minimum). Since both ef- 
fective energy differences, AE^%j^j- (Eq. ([32])) and effec- 
tive potential UeS (Eq. (f22|) ). govern the diffusion of the 
ad-atom, and since t/eff quadratically depends on the am- 
plitude of the force F, we can reasonably expect the same 
quadratic dependence for AE^^yj. However, note that 
while [/off does not take into account the stochastic force 
£_{t) nor the effective crystalline potential, AEg^y^ im- 
plicitly takes into account all these contributions. Figl9] 
exhibits the dependence of AE'g^^ as a function of F^ for 
7 = 2, ^0 = 3.0 and for D = 0.05, 0.1 and 0.2. AE^%jv 
has been calculated from P{x) curves similar to the ones 
of Fig. [H We did not succeed to calculate Aii^g^-^^ as a 
function of F for D < 0.05: the average time needed by 
the particle to escape from an effective crystalline poten- 
tial being too long to obtain good statistics, due to the 
activated character of this event with a 1/ D exponential 
dependence. The linear dependence of Aii^g^^ as a func- 
tion of F^ is observed in the small forces regions where 
the perturbative calculations of Sect. Ill Al scientifically 
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FIG. 9. Effective energy difference Ai^sAw ^^ ^ function of 
ratio ipo/F calculated from histograms Fix). Histogram P(x) 
are calculated from 100 trajectories of 200000 time units for 
each values of ipo/F by solving Eq. (|3H) with 7 = 2.0, tpo = 
3.0, £> = 0.001, 0.01, 0.1. 



sound. We note that increasing the diffusion coefficient 
or the friction 7 (not shown) yields a smaller St AW struc- 
turing effect i.e. as already mentioned, the thermal noise 
challenges the St AW structuring effect. 

The results of this section demonstrates that the St AW 
amplitude and the temperature, acting on the friction 7 
and the diffusion coefficient D allow to tune the struc- 
turation induced by the StAW. These information are 
essential to identify the key parameters if one wishes to 
experimentally implement the dynamic substrate struc- 
turing effect described here. 



VII. CONCLUSION 

In this work, we have studied the solutions of the mo- 
tion equation of an ad-atom diffusing on a substrate sub- 
mitted to a StAW. To this aim, we have applied some 
standard tools of the non-linear physics to a simplified 
version of the motion equation keeping the most relevant 
terms. Noticeably, an effective potential governing the 
slow dynamics of the ad-atom has been derived. We have 



shown that this effective potential controls the distribu- 
tion of the ad-atom x-position even when considering the 
other additional terms (noise, retarded effects, effective 
crystalline potential) in the ad-atom motion equation. 
We underline the relevance of our model when compar- 
ing calculated distribution probabilities to positions his- 
tograms collected from molecular dynamics simulations^. 

Our present study opens different perspectives. 
In a recent publicatior*^, we have discussed about the pos- 
sibilities to experimentally implement the dynamic sub- 
strate structuring effect. We have noticeably mentioned 
the impossibility to produce StAW with a wavelength in 
the nanometer range on the substrate surface using the 
current available experimental setups: to our knowledge, 
and in the state of the art, the production of surface 
acoustic waves with wavelengths in the 100 nm range is 
possible using optical excitationsSi. However, the dy- 
namic structurating effect does not directly depend on 
the wavelength. Actually, the dynamic effect is expected 
to exist as long as the ad-atom does not have the time to 
diffuse on a wavelength A during a period of the StAW 
i.e. 



A;< 



D, 



eff 



(33) 



where Deg is the effective diffusion coefficient of the ad- 
atom in the effective crystalline potential. One may 
thus consider the possibility to use StAW with hundred 
nanometers wavelengths. Our present model and study 
will then be a very fruitful tool to evince the optimized 
parameters (especially the temperature and the StAW 
amplitude) leading to an efficient structuring effect. 

Theoretically, Eq. ([TT]) exhibits a very rich and com- 
plex dynamics. We have studied the solutions of this 
equation, but we currently consider the possibility to de- 
rive the Fokker-Planck equations for the ad-atom prob- 
ability density functions. Such equations would allow to 
directly derive the position distribution P{x), and even 
perhaps to derive the effective diffusion coefficient de- 
scribing the diffusion of the ad-atom between the differ- 
ent basins of the effective potential. 

Finally, in this paper, we have focus on the structuring 
effect induced by the StAW on the position distribution 
of the ad-atom and have eluded the study of the dynam- 
ics of the ad-atom. We have nevertheless exhibited the 
possibility for the ad-atom to follow different qualitative 
types of trajectories depending on the parameters 7, F 
and (/3o- The study of the ad-atom dynamics, including 
the dependence of the StAW force on the substrate lattice 
parameter scale is a natural perspective to this work. 
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